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ABSTRACT 

We discuss the implications of a stellar-dominated UV background at high redshifts for the star forma- 
tion history of Lyman-break galaxies (LBGs) and the thermal and ionization state of the intergalactic 
medium (IGM). The composite spectrum of 29 LBGs evaluated by Steidel et al. (2000) at (z) = 3.4 
can be well fit by a stellar population with ongoing star formation, a Salpeter initial mass function, 
modest or negligible dust reddening, and no intrinsic Hi photoelectric absorption. Fading starbursts 
in which star formation has ceased for 10 7 yr or more cannot reproduce the observed flux shortward 
of 1 Ryd. The small Hi optical depth in LBGs suggests that the neutral gas from which stars form is 
most likely contained in compact clouds of neutral gas with small covering factor. The escape fraction 
of H-ionizing photons must be close to 100 percent for the observed sample of LBGs. The spectrum of 
ionizing photons produced by a stellar population with ongoing star formation is similar to that of QSOs 
between 1 and 3 Ryd, but becomes softer between 3 and 4 Ryd and drops sharply shortward of 4 Ryd. 
A galaxy-dominated UV background appears inconsistent with the observed Hen /Hi opacity ratio at 
z = 2.4, but might be able to explain the Siiv/Civ abundances measured at z > 3 in quasar absorption 
spectra. A scenario may be emerging where star-forming galaxies reionized intergalactic hydrogen at 
z > 6 and dominate the 1 Ryd metagalactic flux at z > 3, with quasi-stellar sources taking over at lower 
redshifts. If the large amplitude of the H-ionizing flux estimated by Steidel et al. is correct, hydrody- 
namical simulations of structure formation in the IGM within the cold dark matter paradigm require a 
baryon density (to explain the observed opacity of the Lyct forest in quasar absorption spectra) which 
is similar to or larger than that favoured by recent CMB experiments, and is inconsistent with standard 
nucleosynthesis values. 

Subject headings: cosmology: theory — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 



['he strength and spectrum of the UV ionizing back- 
ground at high redshift and the nature of the sources 
j^iich reionized the hydrogen component of the intergalac- 
pe medium (IGM) are two of the big outstanding questions 
! tITj observational cosmology. QSOs and star- forming galax- 
rfei have long been considered the two prime candidates, 
^jhas been argued that QSOs may fall short of producing 
a metagalactic flux as high as that inferred from the prox- 
imity effect at high redshift (Bechtold 1994; Giallongo et 
al. 1996; Cooke et al. 1997; Scott et al. 2000) since their 
space density declines rapidly at early epochs (Shapiro & 
Giroux 1987; Miralda-Escude & Ostriker 1990; Meiksin 
& Madau 1993; Madau et al. 1999). The detection of a 
numerous population of Lyman-break galaxies (LBGs) at 
z w 3 (Steidel et al. 1996) makes the idea of massive stars 
dominating the UV background at early epochs quite plau- 
sible. At 1500 A their emissivity outweighs that of QSOs 
by about a factor of fifteen (Madau et al. 1999). Until 
recently it seemed likely that only a small fraction of the 
ionizing photons emitted within starburst galaxies could 
escape the dense Hi layers into the IGM, as is observed in 



their low-redshift counterparts (Leitherer et al. 1995; Hur- 
witz et al. 1997). The detection by Steidel et al. (2000) 
of flux beyond the Lyman limit (with no significant indi- 
cation of a break at 912 A) in a composite spectrum of 29 
LBGs at (z) = 3.4 thus comes as a surprise. 

In this Letter we assess the significance of the Stei- 
del et al. (2000) findings by comparing their spectrum 
to model spectra calculated with the population synthesis 
code STARBURST99 (Leitherer et al. 1999). We discuss 
the implications of this result for the star formation his- 
tory and gas content of these galaxies, for the spectrum 
of the UV background, for the cosmological reionization, 
and for the baryon content of the Universe. 

2. IMPLICATIONS FOR THE STAR FORMATION HISTORY 



The most striking feature of the composite spectrum of 
Steidel et al. is the apparent lack of a significant drop 
at the Lyman limit. Stars with effective temperatures 
< 30, 000 K have a strong intrinsic break at the Lyman 
edge. The composite spectrum shortwards of 1 Ryd is thus 
dominated by O stars with masses J> 20M Q and lifetimes 
<J 10 7 yr. Note that the signatures of O stars have already 
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been detected at longer wavelengths (Pettini et al. 2000) so 
their presence does not come as a surprise. Considering the 
large space density and luminosities of LBGs (Steidel et al. 
1996) it seems very unlikely that these galaxies have been 
caught during a very short starburst. In Fig. 1 we compare 
the observed spectrum with a model spectrum calculated 
with the population synthesis code STARBURST99 (Lei- 
therer et al. 1999). We have assumed a constant rate of 
star formation lasting for t* = 10 8 yr, a Salpeter initial 
mass function (IMF) from 1 to 100 M Q , and solar metal- 
licity. The spectrum is dominated by stars in the mass 
range 15 — 60 M Q , and has been convolved with the model 
of the opacity distribution of the Lya forest clouds and 
Lyman-limit systems of Madau (1995). To improve the fit 
to the data a modest amount of reddening was introduced, 
Eb-v = 0.03 for a SMC (absorption only, no scattering) 
extinction curve (Pei 1992). 

The agreement between the model spectrum and the 
data is rather good, although the observed flux shortward 
of 1 Ryd appears somewhat higher than in the model. 
The discrepancy is, however, within the uncertainties in 
the spectral synthesis and intergalacitc opacity modeling. 
For comparison, we also show the spectrum of an instan- 
taneous starburst after 10 7 yr. This has basically no flux 
shortward of 1 Ryd. Fading starbursts with short duty 
cycles may thus be problematic (e.g. Kolatt et al. 1999). 
Note that there is no hint of a significant H i photoelectric 
optical depth in the local interstellar medium (ISM). 
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Fig. 1. — Thin solid curve: composite spectrum of 29 LBGs at 
(z) = 3.4 from Steidel et al. (2000). Thick solid curve: synthetic 
spectrum of a stellar population with continuous star formation and 
Salpeter IMF. A modest amoung of reddening with a SMC extinc- 
tion curve and Eb—v = 0.03 was included. Dashed curve: synthetic 
spectrum of an instantaneous starburst after 10 7 yr. The model spec- 
tra were calculated with STARBURST99 (Leitherer et al. 1999). 

3. THE NEUTRAL HYDROGEN DISTRIBUTION AND 
ESCAPE FRACTION OF LYMAN-BREAK GALAXIES. 

If star formation in Lyman-break galaxies is indeed 
continuous and lasts for about a hundred Myr there 
should be a substantial amount of neutral hydrogen 
present in these systems. The LBGs in Steidel et 
al. sample have typical star formation rates of M* = 
10 - 50M Q yr _1 . They should thus contain M « 
10 9 - 5 (30M Q yr-VM*) (i*/10 8 yr) M of cold gas if the 
cold gas to stellar mass ratio were of order unity. Tak- 
ing a typical half-light radius of 2kpc as characteris- 
tic size and assuming that the cold gas is distributed 
homogeneously, the typical density and column density 
should be ~ 3 (M/10 9 5 M Q ) (r/2 kpc)" 3 cm" 3 and 2 x 
lO 22 (M/lO 9 - 5 M )(r/2kpc)- 2 cur 2 , respectively. These 



values are higher than those of damped Lya systems and 
a large Hi optical depth would therefore be expected. 

A luminosity in ionizing photons of at least ~ 
lO 44 - 5 (M/lO 9 5 M )(r/2kpc)" 3 ergss- 1 is needed to keep 
this large amount of gas photoionized. The required lumi- 
nosity is a factor 10 — 30 higher than that of the LBGs in 
the sample. Lyman-break galaxies should therefore either 
have a small cold/neutral gas mass to stellar mass ratio in 
their luminous regions or the cold gas must be contained 
in small compact regions with small covering factor. 



4. THE UV ESCAPE FRACTION OF GALAXIES AND THE 
REDSHIFT EVOLUTION OF THE UV BACKGROUND 

The absence of a significant intrinsic H i photoelectric 
opacity in the spectra of LBGs implies that those ionizing 
photons which avoid absorption by dust can escape the 
galaxy unimpeded. After correction for intergalactic ab- 
sorption, the mean ratio of emergent flux density at (rest- 
frame) 1500 A to 900 A is equal to 4.6±1 according to Stei- 
del et al. (2000). Assuming this ratio is typical of all LBGs 
at z w 3 (this might not be true since the composite spec- 
trum is drawn from the bluest quartile of the LBGs popu- 
lation) Steidel et al. estimate J21 = 1.2 ±0.3 for the galaxy 
contribution to the ionizing flux at the Lyman limit (here 
J21 is measured in units of 10~ 21 ergss -1 cm~ 2 Hz 1 sr _1 ). 
This is about a factor four larger than the contribution of 
QSOs at this redshift, and may solve the problem of the 
missing sources for H-reionization (Madau et al. 1999). 
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Fig. 2. — The Hell line blanketing opacity as a function of red- 
shift for a UV background dominated by QSOs and stars. The dif- 
ferent curves show the optical depth for a stellar contribution to 
the flux at 1 Ryd equal to (0, 0.5, 1.5, 2.5 and 5) times that of 
QSOs (from bottom to top). The Hi optical depth is fixed to 
Thi = 0.0038(1 + z) 3 - 46 , while the ionizing flux at z = 2.4 from 
QSOs only is J21 = 0.3 in all cases (Haardt & Madau 1996). The 
curves correspond to an escape fraction of Lyman-continuum pho- 
tons from star-forming galaxies of (0, 10, 30, 50 and 100) percent. 
The data point is from Davidsen et al. (1996). 

Could the UV background be dominated by stars at all 
redshifts ? Comparison of both H 1 and He 11 absorption 
in quasar spectra can be used to constrain the relative 
contribution of massive stars and QSOs to the metagalac- 
tic flux. While in three out of four of the observed QSO 
spectra with a He 11 Lya forest the interpretation is made 
difficult by the large fluctuations measured in the Hen 
scattering opacity - perhaps the signature of a delayed 
Hen — > He in rcionization at z > 3 - (see, e.g., Reimers 
et al. 1997; Heap et al. 2000) the spectrum of HS1700+64 
at somewhat lower redshifts (2.4 < z < 2.7; Davidsen et 
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al. 1996) shows smaller fluctuations and allows a more re- 
liable determination of the average Hen optical depth, 
t'hcII = 1.0 ± 0.07. This value is in good agreement 
with that expected for a quasar-dominated ionizing back- 
ground filtered by the IGM (Haardt & Madau 1996) with 
a frequency dependence of the intrinsic quasar spectrum, 
/„ oc i/- l s (Machacek et al. 2000; Theuns et al. 1998; 
Croft et al. 1997). In Fig. 2 we show the Hen line blan- 
keting optical depth for a UV background with different 
ratios of QSO and galaxy contribution, for a given Hi 
scattering opacity. The five curves correspond to an es- 
cape fraction of H-ionizing photons from galaxies of (0, 
10, 30, 50, 100) percent (bottom to top). For fixed Hi 
scattering opacity the Hen optical depth increases with 
increasing stellar contribution to the total flux. 

A UV background dominated by stars is inconsistent 
with the observed ratio of He n and H i optical depth at 
z ~ 2.4. The average escape fraction of LBGs at this red- 
shift must be <^ 10% similar to the small / esc _ between 
5 and 10% - inferred in local starbursts (Leitherer et al. 
1995; Hurwitz et al. 1997) but very different from the close 
to 100% of the Steidel et al subsample at z ~ 3.4. If the 
results of Steidel et al. were representative of the entire 
LBG population, the average escape fraction would have 
to decrease very rapidly with decreasing redshift. 
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Fig. 3. — Extrapolated far-UV synthetic spectrum of a stellar pop- 
ulation with continuous star formation. A SMC reddening curve as 
modeled by Pei (1992) and the dust opacity distribution of LBGs as 
given by Steidel et al. (1999) have been assumed. The line indicates 
the slope of a typical QSO spectrum with /„ oc y -18 . 



5. A UV BACKGROUND DOMINATED BY STELLAR 
SOURCES AT HIGH REDSHIFT 

Having found a simple model which reproduces the ob- 
served spectrum in the wavelength range 870 A to 1600 A 
we can now extrapolate the UV spectrum to shorter wave- 
length. The result is shown in Fig. 3. We have assumed 
the full range of dust optical depth as reported by Steidel 
et al. (1999) with no dependence of dust optical depth on 
galaxy luminosity, together with a SMC reddening curve 
(absorption only, no scattering; Pei 1992). We also as- 
sumed that the null detection of Hi absorption due to 
the ISM in the Steidel et al. sample is representative of 
all LBGs. Just beyond the Lyman limit the spectrum is 
very similar to that of a typical QSO, somewhat harder 
than the canonical —1.8 power law (Zheng et al. 1998). It 
is only at wavelengths shorter than 400 A that it softens 
significantly compared to a quasar spectrum. Beyond the 
He ii Lyman edge at 4 Ryd there is a strong break. 

The similarity to a QSO spectrum just beyond 1 Ryd 



may explain the already high temperatures of 10 4 K at 
z > 3 (Haehnelt & Steinmetz 1998; Schaye et al. 2000). 
A UV background which is progressively dominated by O 
stars may also be the solution to the increasing Siiv/Civ 
ratio with increasing redshift observed in metal absorp- 
tion systems of intermediate column densities (Songaila & 
Cowie 1996; Giroux & Shull 1997; Boksenberg et al. 1998). 

How different could the integrated stellar spectrum be? 
The upper cut-off of the IMF and the metallicity will 
strongly affect how far the spectrum extends into the far 
UV (Tumlinson & Shull 2000; Bromm et al. 2000; Oh et 
al. 2000). We have assumed an upper cut-off of 100M Q 
and solar metallicity. The dust opacity in the far UV de- 
pends on smaller grains than usually probed by reddening 
analysises and is rather uncertain. In the the model of Pei 
et al. (1992) for the SMC opacity used here, the dust opac- 
ity drops shortwards of 1000 A. This makes the spectrum 
bluer not redder at wavelengths shorter than 1000 A. The 
stellar spectrum of massive stars with strong mass loss are 
complex and require non-LTE calculations including the 
effect of stellar wind outflow and spherical extension. In 
particular, the strength of the 4 Ryd break depends sig- 
nificantly on metallicity, mass loss rates and details of the 
modelling (Kudritzki 2000; Pauldrach et al. 2000). 

6. THE BARYON FRACTION OF THE UNIVERSE 

One interesting consequence of the large amplitude of 
the UV background suggested by Steidel et al. is the im- 
plied conflict of the Lya opacity measurements in QSO 
spectra with the nucleosynthesis constraint for the baryon 
density. The optical depth for Lya scattering of the pho- 
toionized IGM scales as (Rauch et al. 1997a) T^ ya oc 
i?(z)- 1 (fi b /i 2 ) 2 7 1 ~°- 7 r~ 1 , where T is the photoionization 
rate, T is the gas temperature, H{z) is the Hubble con- 
stant at redshift z, Qb is the baryon density parame- 
ter, and h is the present-day Hubble constant in units 
of 100 kms -1 Mpc~ x . The investigation of Rauch et al. 
implies that for a nucleosynthetic value of fl^h 2 = 0.019 
(Buries & Tytler 1998) the Lya opacity can be matched if 
r ~ 5 x 10 s . For a spectrum f u oc v~ 1B this corre- 
sponds to an ionizing flux J21 = 0.12. The large value of 
J21 = 1.2 ± 0.3 suggested by Steidel et al. implies instead 
ilbh 2 ~ 0.06. The D/H measurements are therefore in 
strong conflict with the Lya opacity measurements if both 
the standard nucleosynthesis calculations and the value 
of J21 suggested by Steidel et al. are correct. This dis- 
crepancy may be even stronger if there were a significant 
contribution to the H-ionizing background from faint Lya 
emitters with weak stellar continuum flux as suggested by 
Kudritzki et al. (2000). 

In Fig. 4 we have compiled from the literature a number 
of values for fl\>h 2 implied by the measured Lya opacity 
assuming J21 — 0.6 as may be appropriate if LBGs with 
redder spectra than those in the Steidel et al. sample had 
significant smaller escape fraction for H-ionizing photons. 
The values at z — 3 are taken from Rauch et al. (1997b), 
Gnedin (1998), Nusser & Haehnelt (2000), Theuns et al. 
(2000), and at z = 2 from Cen et al. (1994) and Hern- 
quist et al. (1996). We have scaled all values to a tem- 
perature of 15000 K at z = 3 and 11000 K at z = 2 as 
suggested by the temperature determination of the IGM 
by Schaye et al. (2000). All values for fib scale oc J 21 and 
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oc T - 35 . We also show the nucleosynthesis constraint from 
the D/H measurement as given by Buries & Tytler (1998) 
and the combined constraints from the CMB experiments 
BOOMERANG and MAXIMA (Jaffe et al. 2000). 
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Fig. 4. — The baryon fraction inferred from the Lya opacity as 
described in the text. The values inferred from D/H measurements 
(Buries & Tytler 1998) and CMB experiments (Jaffe et al. 2000) are 
also shown. 



7. DISCUSSION AND CONCLUSIONS 

We have shown that the detection of flux beyond the 
Lyman limit from LBGs is most easily explained as being 
due to the O stars present in a star forming region with 
ongoing star formation, a Salpeter IMF, modest or negligi- 
ble reddening, and negligible absorption due to Hi in the 
local ISM. Fading starbursts in which star formation has 
ceased for 10 7 yr or more cannot reproduce the observed 
flux shortward of 1 Ryd. The small H I optical depth can- 
not be due to photoionization. The cold neutral gas from 



which they stars form must be contained in compact small 
regions with small covering factor. 

The resulting ionizing background will be dominated by 
O stars shortward of the Lyman limit and is as hard as a 
QSO spectrum between 1 and 3 Ryd. Shortward of 3 Ryd 
the spectrum is softer than a QSO spectrum and there 
is strong drop shortward of 4 Ryd. A consistent picture 
seems to emerge where the hydrogen of the Universe is 
reionized by stars at z > 6 while the reionization of helium 
is due to QSOs and is delayed to z ~ 4 or 3. The contri- 
bution of stars to the ionizing UV background is equal or 
larger than that QSOs at least up to z ~ 3. At smaller 
redshift the ratio of Hen and Hi opacity and the ob- 
served Siiv/Civ ratios in the associated metal absorption 
systems of intermediate column density QSO absorption 
systems argue for a UV background dominated by QSOs. 
This requires a substantial decrease of the escape fraction 
of ionizing photons from star-forming galaxies with de- 
creasing redshift. The amplitude of the UV background in- 
ferred from the composite spectrum of Steidel et al. makes 
the measurements of the opacity of the Lya forest in QSO 
absorption spectra inconsistent with D/H measurements 
if the spectrum is representative for all LBGs and nucle- 
osynthesis calculations for D/H are correct. The inferred 
value for fib is similar to that that indicated by the CMB 
experiments Boomerang and Maxima if J21 = 0.4 — 0.6. 
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